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We investigate the localization properties of atoms moving in a three-dimensional optical lattice 
in the presence of an uncorrelated disorder potential having the same probability distribution P{V) 
as laser speckles. We find that the disorder-averaged (single-particle) Green’s function, calculated 
via the coherent potential approximation, is in very good agreement with exact numerics. Using the 
transfer matrix method, we compute the phase diagram in the energy-disorder plane and show that 
its peculiar shape can be understood from the self-consistent theory of localization. In particular, 
we recover the large asymmetry in the position of the mobility edge for blue and red speckles, which 
was recently observed numerically for correlated speckle potentials. 
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I. INTRODUCTION AND MODEL 
HAMILTONIAN 


Recently there has been a growth of experimental and 
theoretical studies [l| on Anderson localization of ultra¬ 
cold atoms exposed to disordered optical potentials, in¬ 
cluding quasi-periodic lattices and blue-detuned laser 
speckles. The latter behave quite differently from typi¬ 
cal models of randomness considered for condensed mat¬ 
ter systems. Indeed, the potential distribution P{V) of 
blue speckles is not Gaussian but follows the exponential 
(Rayleigh) law ii: 


P{V) 


0(U + Uo) 



( 1 ) 


where 0 is the Heaviside function and Vq is the disorder 
strength. In Eq. (U), we have shifted the potential by its 
average value so that V = 0. This distribution is bounded 
from below by —Vq and is asymmetric, implying that odd 
moments of the potential will be non zero. The proba¬ 
bility distribution of red-detuned potentials, which have 
not yet been implemented experimentally, is obtained by 
simply changing V to — U in Eq. O- A second funda¬ 
mental aspect of optical speckles is that they are spatially 
correlated, with a typical grain size of the order of 1 /rm. 

Three different experiments [3-@| claimed the obser¬ 
vation of 3D Anderson localization of atoms exposed to 
blue speckles. However, the estimated mobility edge - 
namely the critical value A = Ac of the energy separat¬ 
ing localized (A < Ac) from extended (A > Ac) states - 
is in all cases larger (and even much larger in @ 
the current theoretical and numerical predictions 
The question arises whether the problem is due to inade¬ 
quate experimental measurements (it is for example very 
difficult to correctly analyze the density profiles after a 
relatively short time expansion in the presence of disor¬ 
der [il) or to incorrect theoretical predictions for the 
mobility edge. 

From the theory side, exact numerical calculations of 
the mobility edge obtained via the transfer-matrix tech¬ 



nique 0 have revealed a discrepancy with previous esti¬ 
mates based on the self-consistent theory of Anderson lo¬ 
calization (SGTL) (MU. In particular it was found that 
the on-site potential distribution is extremely important. 
Indeed there is a large difference between blue-detuned 
speckles, where the mobility edge is systematically neg¬ 
ative (i.e. below the average potential), and red-detuned 
speckles, where it changes from negative at low Vb to 
positive at large Vq Q- On the other hand, the details 
of the spatial correlation function - beyond the correla¬ 
tion length which fixes the characteristic quantum energy 
scale - do not play a major role. The aforementioned dis¬ 
crepancy with SGTL predictions originates most proba¬ 
bly from the lack of accuracy in the calculation of the 
disorder-averaged single-particle Green’s function, which 
is a fundamental ing redient of the SGTL approach. In¬ 
deed, in Refs. lol-Hll this quantity has been computed 
within the self-consistent Born approximation (SGBA), 
which is by construction insensitive to the color of the 
speckle, and can only apply for a sufficiently weak disor¬ 
der. 

In this work we disentangle the role of the on-site po¬ 
tential distribution from the effects of spatial correla¬ 
tions by studying a spatially uncorrelated 3D Anderson 
model with a Rayleigh potential distribution. We cal¬ 
culate the Green’s function numerically, and show that 
the obtained results are well reproduced by the coherent- 
potential approximation (GPA), which takes into account 
the full statistical properties of the potential. The GPA 
self-energy is then injected in the SGTL to estimate the 
position of the mobility edge, which turns out to be 
in reasonably good agreement with the exact transfer- 
matrix calculations. It is worth mentioning that the GPA 
method has not yet been generalized to spatially corre¬ 
lated speckle potentials, although this generalization has 
been attempted for other types of spatially-correlated 
disorder 

In first quantization, the 3D Anderson model is given 
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by 

H = '£-J\i){j\ + '£v,\t){il (2) 

{ij) » 

where indices i,j label the sites of the lattice, J is the 
hopping term between nearest neighbors, and Vi the ran¬ 
dom external potential which is assumed to be uncorre¬ 
lated, ViVj = {V^)5ij, and obey the potential distribution 
in Eq. O- For convenience, in the following we will use 
units J = 1 and h= 1. For the specific case of the Ander¬ 
son model, there is no need to study separately the blue 
and red-detuned speckles. Indeed, the cubic lattice with 
nearest-neighbor coupling being bipartite, the sign of J is 
irrelevant and does not affect localization properties. As 
a consequence, a change V —> — F is equivalent to revers¬ 
ing the sign of the Hamiltonian, H —>■ —H. This means 
that all results obtained here for the blue speckle with 
the Rayleigh distribution, including the phase diagram, 
apply also to the red speckle under the change E ^ —E, 
as shown in Fig. [2] In other words, by studying the two 
mobility edges around E = —6 and E = -|-6, one covers 
both the blue and red-detuned cases respectively. 

The plan of the paper is as follows. In Sec. m we 
present the exact phase diagram for the 3D Anderson 
model with the Rayleigh potential distribution, Eq. o, 
obtained numerically via the transfer-matrix technique. 
In Sec. imi we calculate the disorder-averaged single¬ 
particle Green’s function numerically, and analytically 
using the CPA and the SCBA. The obtained results are 
then used in Sec. m to estimate the position of the mo¬ 
bility edge via the SCTL and to compare it with the 
exact numerical results of Sec. HD Finally, we examine 
in Sec. E how our findings help to interpret numerical 
and experimental results for cold atoms in a continuous 
speckle potential. 


II. EXACT PHASE DIAGRAM 

The transfer-matrix method consists in computing re¬ 
cursively the conductance of a bar with transverse sec¬ 
tion M X M sites, and length L ^ M. The logarithm of 
the conductance (averaged over disorder realizations) de¬ 
cays linearly with the length of the bar, giving access to 
the quasi-ID localization length Am- The scaling theory 
of localization indicates that Xm/M decreases with M 
in the localized regime, increases in the diffusive regime, 
and is constant for large M at the mobility edge [l^ . The 
various Xm/M vs. energy curves - for various M values 
- thus cross at the critical energy, the mobility edge, as 
shown in Fig.[I] By gathering results at various values of 
the disorder strength Vq, we obtain the phase diagram of 
the Anderson model, shown in Fig. [2] for a blue-detuned 
speckle. 

We note that the trajectory of the mobility edge be¬ 
haves very differently at the two band edges. Indeed, 
starting from the left band edge of the clean model, 



E 


FIG. 1. (Color online) Transfer-matrix results for the 3D 
Anderson model with uncorrelated, blue-detuned speckle po¬ 
tential of strength Vo = 8. Each curve displays the ratio of the 
localization length Am of a long bar with cross-section M x M 
to the bar transverse size M, as a function of energy. The 
various curves from M = 16 (least steep curve) to M = 31 
(steepest curve), cross at the mobility edge Ec ~ —9.93. 


E = —6, the mobility edge shifts to lower and lower en¬ 
ergy as Vo increases until Vq ~ 10, which is close to the 
critical disorder needed to localize all states. Since the 
potential distribution has a sharp cut-off at V = — Vq, 
there are strictly no states below E = —6 — Vo- In con¬ 
trast, the spectrum is not bounded from above. At the 
right band edge, E = -|-6, however, the mobility edge 
moves only slightly to the right for very weak disorder, 
reaching if = 6.15atVo = l, and then moves backward 
in a monotonous way. This means that even for mod¬ 
erate disorder, starting from Vb ~ li the asymmetry of 
the Rayleigh probability distribution P{V) becomes im¬ 
portant and must be taken into account. This rules out 
the use of the SCBA for the calculation of the disorder- 
averaged Green’s function, as we will show below. 


III. DISORDER-AVERAGED GREEN’S 
FUNCTION 

In this section we calculate the disorder-averaged 
single-particle Green’s function numerically, and analyti¬ 
cally using two different approximation schemes, namely 
the SCBA and the CPA. We also discuss two related 
quantities, the (disorder-averaged) density of states and 
the effective band edge (neglecting Lifshitz tails). 

In the absence of disorder {Vo = 0), the Hamiltonian 
E[ becomes diagonal in momentum space: E[ = EIq = 
e(k)|k)(k|, where e(k) = —2(cosfca; + cosky + cos kz) 
is the energy dispersion of the atom in the cubic lat¬ 
tice. Hence the diagonal part of the free-particle Green’s 
function in configuration space representation Go{E) = 
{n\{E — i/o)~^|n) is translationally invariant and given 
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FIG. 2. (Color online) Exact localization phase diagram for 
the 3D Anderson model with blue-detuned speckle potential 
obeying the Rayleigh distribution, Eq. m, as obtained from 
the transfer-matrix technique and finite-size scaling. Notice 
that there are rigorously no states below the solid line corre¬ 
sponding to E = —6 —Vq. The phase diagram for red-detuned 
speckles is simply obtained under the change E —>■ —E, as 
shown in the inset. 


by 


GoiE) 


(27r)3 E - e(k) -b iO’ 


( 3 ) 


Following Joyce [l^, we can express the unperturbed 
Green’s function of a cubic lattice as Go{E) = 
P{6/E)/E, where 


P{z) 


1-9^4 

( 1 - 0 ^( 1 + 30 



( 4 ) 


Here ^ and ki are functions of z defined as 

1/2 


00 = 
ki{zf = 


1 - Ql - zV9 

1 + Vl - 
160 


( 5 ) 


( 6 ) 


( 1-001 + 30 ’ 

and K is the complete elliptic integral of the first kind. 


A. Numerical computation and analytical 
estimates of the self-energy 

In the presence of a random potential, the transla¬ 
tional symmetry is restored only after averaging over dif¬ 
ferent realizations. The corresponding disorder-averaged 
Green’s function is given by 

^ E-e(k)-E{E,k) ’ 

where E = E'-biS" is the self-energy. The real part of the 
self-energy E'(£', k) represents how the energy of a plane 


wave with wavevector k is shifted from e(k) (“renormal¬ 
ized”) under the influence of the disorder, while the imag¬ 
inary part E"(FJ, k) - which is always ne gati ve - yields 
the broadening of the energy distribution [17|. 

Numerically computing the self-energy is possible, but 
not completely straightforward. The starting point is the 
temporal representation of the Green’s function as: 


1 

E-H+ i0 



^iEt 


di. 


( 8 ) 


As mentioned previously, the average Green’s function 
in presence of disorder is invariant by translation, i.e. 
diagonal in momentum space. Thus, one has: 


(k|G(i?)|k) = 


E - e(k) - E(F;, k) 

poo _ 

—i / (kje^'-^^lk) e‘'®* dt. 


( 9 ) 


The numerical calculation then amounts to propagating 
an initial plane wave |k) with the disordered Hamiltonian 
E[ during time t (with periodic boundary conditions) and 
to computing the overlap of the time-evolved state with 
(k|, followed by a Fourier transform from time to energy; 
subtracting E — e(k) from the inverse of the result yields 
the self-energy. This procedure is then repeated for sev¬ 
eral independent realizations of the disorder to perform 
disorder averaging. In order to obtain small statistical 
fluctuations, a rather large number (of the order of 10®) 
of disorder realizations is needed. Moreover, a sufficiently 
large system - much larger than the mean free path ~ has 
to be used to reduce hnite-size effects. In our numerical 
results, we found that the self-energy has a very weak de¬ 
pendence on momentum, which is expected for spatially 
uncorrelated potentials, ViVj = VgSij. In Fig. [3l we plot 
our numerical results for the real (upper panel) and imag¬ 
inary (lower panel) parts of the self-energy, calculated at 
k = 0 for two different values of the disorder strength. 

In the following we compare our exact numerics with 
several approximate methods to calculate the self-energy. 
By construction, these methods yield estimates which are 
independent of momentum, E(£i,k) = E(E). The sim¬ 
plest one is the Born approximation, which is given by: 


E{E) = Vo^GoiE). (10) 


A slight, but simple, improvement is the self-consistent 
Born approximation (SCBA) where the disorder-free 
Green’s function in Eq. (HOD is modified self-consistently, 
leading to the equation: 

E{E) = Vo^Go{E-E). ( 11 ) 


This quantity can be easily calculated numerically by suc¬ 
cessive iterations, starting from the Born approximation, 
Eq. (ITU)) . Since only the second moment of the potential 
distribution appears in the rhs of Eq. (fTTT) . the SCBA 
self-energy does not depend on the details of the poten¬ 
tial distribution P{V), which makes sense only for suffi¬ 
ciently weak disorder. 






















A better approximation scheme at stronger disorder 
is given by the CPA [l^. The basic idea of the CPA 
is as follows (for a pedagogical discussion see Ref. [H): 
one isolates a single site i where the potential is chosen 
randomly according to P{V), and replaces the surround¬ 
ing sites by an effective homogeneous medium character¬ 
ized by a uniform self-energy. One can then compute the 
single-site t-matrix of site i embedded in the surrounding 
effective medium. The consistency condition is that the 
t-matrix, averaged over the potential distribution at site 
i, vanishes; this translates in the self-consistent equation 


/ 1 - (P - ^{E))gJ^E - E{E)) 

which clearly depends on the full probability distribution 
of the potential. For weak disorder, the CPA reduces to 
the SCBA, as is clear from Eq. m using the Taylor ex¬ 
pansion 1/(1 — cc) ~ 1 -I- x, with X = {V — T,)Go{E — E). 
As the disorder strength increases, however, higher order 
terms in the expansion become important and all mo¬ 
ments of the potential distribution (cubic, quartic, etc.) 
will start to play a role. Hence, in general, the CPA 
yields more accurate results than the SCBA, as previ¬ 
ously known for the Anderson model with random-box 
potential [2^ [2l| . By substituting the rhs of Eq. ([T]) in 
Eq. (I12p and performing the integration over P, we ob¬ 
tain 


e '^{Ei{u)-in) = gVo, (13) 

where g = Go{E - E), m = (1 -b ^(Po -b E))/(gPo) 
and Ei(u) is the exponential integral function defined as 
Ei{u) = —PV dte~*/t, PV being the principal value. 
Equation (1131) can be easily solved numerically by a root¬ 
searching algorithm for the complex variable E or, again, 
by an iteration scheme. 

The SCBA and CPA results for the self-energy are 
shown in Eig. [3] as dotted and solid lines, respectively. 
The most striking observation is that the CPA gives 
excellent predictions, almost in perfect agreement with 
our numerical results. Hence, this approximation can be 
safely used in the SCTL to compute the phase diagram 
of the Anderson model with a potential distribution typ¬ 
ical of optical speckles. In contrast, the SCBA is a rather 
poor approximation, even at moderate disorder strength. 

By construction, in the SCBA the real (imaginary) 
part of the self-energy is an odd (even) function of the en¬ 
ergy and of the potential strength. We see in Eig. |3] that 
this feature disappears in the CPA, due to the asymme¬ 
try of the Rayleigh distribution, Eq. dH). In particular, 
for a fixed disorder strength, the imaginary part of the 
self-energy is much larger (in modulus) around E = 6 a.s 
compared to E = —6, suggesting that disorder scattering 
is much stronger in red speckles. 

The dependence of the self-energy on the speckle color 
is most evident in the atomic limit, either for Vq ^ I 
or A ^ I, where all states are deeply localized. By 
neglecting the tunneling term in the Hamiltonian, one 



FIG. 3. Energy dependence of the real (E') and imaginary 
(E”) parts of the self-energy calculated using the coherent 
potential approximation (solid line, Eq. (I12ll l and the self- 
consisteirt Borir approximatioir (dotted line, Eq. (1111) 1 for two 
values of the disorder strength: Vb = 1 (left) and Vb = 5 
(right). Also shown is the atomic limit approximation (dashed 
line, from Eqs. usiiini)) for Vb = 5. In both cases, the CPA is 
a truly excellent approximation to the exact numerical results 
(for k = 0) plotted as circles, reproducing almost all details 
of both the real and imaginary parts of the self-energy (note 
that for the numerical computations, system sizes of with 
M — 28 were used). In contrast, the SCBA is a rather poor 
approximation, even at moderate disorder strength Vb = 1 
where it is supposed to work best. 


is left with a single-site problem. The disorder-averaged 
Green’s function then takes the simple form 

By substituting the rhs of Eq. dl]) in Eq. ((Ill) and per¬ 
forming the integration over V, we find 


G(E) 




-ine-^/^°0(E + Vo) , 


(15) 


where 


r+°° p -2 

f(x) = PV / -dz. (16) 

J-i x-z 

The self-energy, calculated from Eqs. dm mi), is shown 
in Eig. dl for Vq = 5 as a dashed line and we see that it 
agrees with the exact numerics at high energies. 

We stress that the exact energy spectrum is bounded 
from below by E = — 6 — Vq, which implies that the imag¬ 
inary part of the self-energy (E") is strictly zero below 
this value. However, numerical results show that there 
are very few states immediately above E = — 6 — Vq, in 
the so-called Lifshitz tail, where E" is extremely small. 
It only raises at a significantly larger energy E k, —8.4 
for Vo = 5 (while -6 - Vo = -H), see Fig. [3] or Fig. 1 
below. Thus, there is an effective band edge which is 

























higher than the rigorous band edge. Both the SCBA and 
the CPA, being insensitive to Lifshitz tails, have a band 
edge below which E" vanishes. We see in Fig. [3]that the 
CPA band edge is in excellent agreement with exact nu¬ 
merics, while the SCBA band edge is displaced towards 
negative energies. 


B. Density of states and band edge 

In Fig. m we plot the averaged density of states (DOS) 
p = —l-mG{E)/'K calculated within the CPA for increas¬ 
ing values of the disorder strength, compared to ex¬ 
act numerical values obtained from the self-energy using 
Eq. ([T]). Again, the agreement between the CPA and ex¬ 
act numerics is very good. As disorder increases, the peak 
in the DOS shifts towards negative energies and becomes 
strongly asymmetric [l^. In particular, the DOS devel¬ 
ops exponential tails (decaying as exp(—A/Vb)) at high 
energies, whereas at negative energies it vanishes more 
and more sharply with increasing Vq. From Eq. dig , we 
indeed see that in the atomic limit 

p(E) = ^e-^/^°Q{E + Vo), (17) 

eVo 

which coincides with the potential distribution, Eq. o, 
after replacing V hy E . As shown in Fig.lU this formula 
correctly reproduces the high-energy tails of the DOS. 
We emphasize that Eq. (HU can be understood directly 
from Eq. (|T4)) using the formula (x -|-i0)“^ = PV(l/a;) — 
i7r^(a:). We mention that the density of states of atoms 
in spatially-correlated speckle potentials has been studied 
both numerically and analytically in one dimension [^ . 

We now turn our attention to the effective band edge 
and derive analytical formulas for the SCBA and the 
CPA predictions for weak disorder. The unperturbed 
Green’s function is first expanded around the unper¬ 
turbed band edge E = —6, yielding Gq{E) = Go(—6 -I- 
E) ~A + By/^, where A = Go(-6) = (14^6-^20^3- 
24V2-36)K^l(2-V3)(V3-V2)]/7r^ Ri -0.252731 and 
B = l/(47r) Ri 0.079577. Substituting this expression in 
Eq. (fTTl) gives the self consistent equation: 

S = (a -f S\/e - a) . (18) 

For weak disorder, we can replace the self-energy in the 
rhs of Eq. (ITOl) by its leading term, E = VqA, which is 
real. As a consequence, the self-energy can only become 
complex if the argument of the square root in Eq. CHD 
becomes negative, that is A = E. This means that the 
SCBA band edge is given by ^ _Q_^AV2_^o{Vf^). 

It is possible to get the next order term - note that by 
construction of SCBA, only even powers of Vq appear - 
using Eq. m- By bringing the term proportional to A 
to the Ihs and taking the square of both sides, we find 
that E satisfies the following quadratic equation: 

E^ - {2AV^ + B2vb‘)E -h A^V^ + B^V^E = 0. (19) 



E 

FIG. 4. (Color online) Density of states calculated in the 
CPA as a function of energy for three different values of the 
disorder strength Vb = 1 (black solid line, first from top at 
A = 0), Fo = 5 (blue) and Vb = 10 (red). The dashed line 
correspond to the strong disorder expression in Eq. (1171) . Dots 
correspond to the exact numerical results for system sizes A'E 
with M = 28 (for Fq = 1), and M = 16 (for Fo = 5,10). The 
CPA is an excellent approximation, reproducing very well the 
existence of an effective band edge at low energy. Only some 
tiny features above the band edge are not correctly repro¬ 
duced. 


The band edge then corresponds to the energy value at 
which the discriminant vanishes. This gives 

£;SCBA ^ _g ^ ^y2 ^ ^ ( 2 O) 

which is shown in Fig. [5] with the red solid line. 

Let us now focus on the CPA band edge. By setting 
-z = gVo, y = gE and x = F/Vq, the self-consistent 
Eq. (fl^ for the CPA self-energy can be written as 



e 

1 — xz + y 


dx. 


( 21 ) 


The argument in the integral is expanded in powers of y 
up to the quadratic term: 



y 

(l-xz)^ 


= fo{z) + fi{z)y + f 2 {z)y^, 


\ 

(1 — xz)^ J 


dx 


( 22 ) 


where fi{z) are functions of 2 defined in the interval 
[—1,0]. For small values of z, f 2 {z) — implying that 
y = F{z) + 0(z®), where 


F{z) 


fojz) 

1 - /l(^)’ 


(23) 


By inserting the asymptotic expansion of the Green’s 
function g ~ A+B\/E — E in the formula E ~ F{gVo)/g 
and Taylor-expanding the rhs, we find with the same level 
of accuracy that 


E(AFo) -b F'(AVo)BVoVe - E 
A + By/E-E 


( 24 ) 
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FIG. 5. (Color online) Band edge Ebe plotted as a function of 
disorder strength, calculated within the SCBA (red squares) 
and the CPA (black circles). Solid lines refer to the corre¬ 
sponding approximate analytical expressions, Eqs. (I2S1) and 
(I25j, obtained for weak disorder. The exact transfer-matrix 
results for the mobility edge are also shown (green diamonds). 


Repeating the same procedure as above, we obtain the 
following approximate formula for the CPA band edge: 


E 


CPA 

BE 


4 — 

'F'{AVo)Vo 

A 


(25) 


F{AVo)Y 

A2 



disorder - between k and E, E = e(k), is no longer cor¬ 
rect in the vicinity of the mobility edge (strong-scattering 
regime), and there is no unique way of defining k{E) and 
i{E). In the simplest approximation, one may assume 
that the dominant effect of strong disorder is to shift the 
band edge and dispersion relation by the real part of the 
self-energy E', so that, for example, k{E) is nothing but 
k[j{E — S'(£l)), where the ‘0’ subscript refers to disorder- 
free quantities. As for the broadening of the spectral 
function (related to the imaginary part E"), it is approx¬ 
imately symmetric and thus expected to have a negligible 
effect on averaged quantities. For lattice models, which 
have an anisotropic dispersion relation, k becomes a vec¬ 
tor quantity so that the loffe-Regel criterion cannot be 
used straightforwardly. A simple generalization has been 
proposed in Ref. [l^, that reads: 

So{E - E’)£{Ef = 47^C'^ (27) 

where Sq{E) is the area of the surface e(k) = E in mo¬ 
mentum space: 

SoiE) = J |Vke|,5(A - e{k))d^k. (28) 

The mean free path can be written as £{E) ^ v{E)t{E), 
where v{E) = Vq{E — E') is the average modulus of the 
particle velocity, defined as 

v^{E) = J |Vkep<5(A - e{k))d^k, (29) 


This formula, shown in Fig. [5] with the black solid line, 
reproduces very well the numerical results (open circles) 
for the CPA band edge at small disorder. By using the 
Taylor expansion E{z) = + 2z^ + 7z'^ + Sdz^ + 0{z^), 

one can obtain the expansion of the CPA band edge in 
powers of Vq: 

F;Cpa = _6 + AV^ + 2A^V^ + + 7A^^ 

-f (2AR2 + 34 ^ 4 ^ + 0(V(f). (26) 

Note that odd-power terms, which are absent in the 
SCBA, appear due to the asymmetry of the Rayleigh 
potential distribution, and that the fourth order term in 
Eq. (1251) is different from the SCBA result, Eq. (1201) . 


IV. SELF-CONSISTENT THEORY OF 
ANDERSON LOCALIZATION 

Starting from the weak-localization corrections to the 
conductivity (or to the diffusion constant), the SCTL 
provides a useful microscopic justification of the loffe- 
Regel criterion ki k. C for the onset of localization in 
3D continuous models, C being a model-dependent con¬ 
stant of order unity. Both the wave number k and mean 
free path i depend on the energy A, but not in any sim¬ 
ple way. Indeed, the usual relation - in the absence of 


and t(A) is the relaxation time due to disorder, which is 
related to the imaginary part of the self-energy as 


1 


Then Eq. (1271) reduces to 


5o(E-E')uo(E-E')2 

E"2 




(30) 


(31) 


Note that for continuous systems with a “massive” dis¬ 
persion relation, E = k^/2m, one has Sq{E) = SirmE 
and Vo{E) = ko/m = ■\j2E/m, so that Eq. m indeed 
reduces to the loffe-Regel criterion. 

Equations dMl [21 ED) allow to extract the complete 
phase diagram of the Anderson model once the CPA 
self-energy has been computed from Eq. (USD, and a 
reasonable choice for the constant C in Eq. has 
been made. For the random-box potential distribution, 
Phox{V) = 0(lF/2 — V)/W, this constant was chosen 
p3 | to reproduce the well-known transfer-matrix result 
W = Wc = 16.5 for the critical disorder strength at the 
center of the band {E = 0). Inserting the box distri¬ 
bution in Eq. (1T21) for the CPA self-energy, one obtains 
Ebox(E = 0) = 0 — i 4.7011 for W = 16.5, which yields, 
using Eq. ®, the constant Cbox = 0.775. For definite¬ 
ness, we will calculate the mobility edge for the uncor¬ 
related speckle potential using the same constant, that 
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FIG. 6. Prediction of the self-consistent theory of localiza¬ 
tion - see Eq. ®1) - for the mobility edge of the Anderson 
model with Rayleigh potential distribution, using the CPA 
self-energy (solid line). The empty circles correspond to the 
transfer-matrix results of Fig. [2l The dashed line represents 
the band edge calculated from CPA, neglecting Lifshitz tails. 


is we set C = Cbox- Other choices of the constant will 
produce qualitatively similar results. 

In Fig.ini we compare the SCTL predictions (solid line) 
with the transfer-matrix results obtained in Sec. m We 
see that the SCTL reproduces the overall shape of the 
exact phase diagram. In particular the generalized loffe- 
Regel criterion, Eq. dan, provides a natural explanation 
for the behavior of the mobility edge at positive energies. 
Indeed we see from Fig. [3] that near E = 6 the imaginary 
part of the self-energy E" is rather large, or equivalently, 
the relaxation time in Eq. (1301) is rather short. As a 
consequence, So{E — S') must therefore increase, which 
means that Ec will move towards the center of the band, 
where So{E) takes its maximum value, 5'o(0) = 92.8756. 

In Fig. [6l one can also see cusps in the SCTL mobil¬ 
ity edge occurring at E = —9.14 for Vq = 11.22 and, 
although it is less evident, at E = —0.83 for Vq = 6.86. 
These features are related to the presence of Van-Hove 
singularities [ 2 ^ in the DOS of the clean system at 
E = ±2. Indeed, at the cusp positions, the real part 
of the self-energy satisfies E — S' = — 2 and E — S' = 2, 
respectively. 

Another interesting feature of Fig. [3] is that the mo¬ 
bility edge at negative energies remains very close to 
the CPA band edge, even for moderate disorder. This 
can be easily understood: In the immediate vicinity of 
the effective band edge E^^^, the density of states and 
dispersion relation are very similar to the ones of the 
disorder-free system near the disorder-free band edge 
E = —6, meaning that they can be obtained from the 
disorder-free quantities by the shift E —?► E — Ebe, where 
Abe = This shift is slightly different from the one 

in Eq. eil),i.e. using S', although the two quantities co¬ 
incide up to order Vq (see above). Using Ebe is a better 
approximation however, as it ensures that the density 


of states exactly vanishes at E = Ebe- Thus, close to 
the effective band edge, one has approximately the aver¬ 
age velocity (the effective mass of the Anderson model is 
m = 1/2): 

vo{E-EBE)^2k{E), (32) 

with k{E) = ^E — Abe, whereas the area of the 
constant-energy shell is 

Sq{E - Ebe) ^ ^-^KEf. (33) 


It also follows that the imaginary part of the self-energy 
in the Born approximation, Eq. OT, behaves like: 


s"(e;) 


^o^fc(A) 

47r 


(34) 


This is in agreement with the square-root behavior of 
E"(E) above the CPA band edge observed in Fig. |31 It 
also means that at low energy the mean free path is pro¬ 
portional to 1/Vq and independent of energy. Insert¬ 
ing Eqs. (I321f34)) in the generalized loffe-Regel criterion, 
Eq. (|3^ . yields 


Therefore, in the weak disorder limit Vq li the effective 
band edge is shifted to the left of A = —6 proportionally 
to Vq , while the distance to the mobility edge is much 
smaller, scaling like Vq. 


V. CONSEQUENCES FOR ANDERSON 
LOCALIZATION OF COLD ATOMS 

A direct implementation of the uncorrelated 3D An¬ 
derson model with cold atoms would require to create a 
tight cubic optical lattice so that only the first band is 
populated, an ultra-cold gas so that ksT is much smaller 
than the bandwidth (proportional to the tunneling rate 
between neighboring sites), and a speckle potential with a 
correlation length much shorter than the lattice spacing. 
Meeting all these requirements in current experiments 
seems rather difficult. 

From the present study, one can nevertheles draw a 
few conclusions about the Anderson localization of mas¬ 
sive particles in realistic optical speckle, as realized in 
recent experiments. Because of finite-range correlations, 
the disorder strength Vq must be compared to the other 
characteristic energy scale, namely the “quantum” cor¬ 
relation energy E^ = a being the correlation 

length. A thorough discussion of the various possible 
regimes can be found in Ref. S Especially, the effective 
band edge and the distance of the mobility edge to the 
effective band edge were found to scale like Vq and Vq re¬ 
spectively, in agreement with our results (with however a 
caveat, see below). In the “quantum” regime, Vq <C E^, 
the de Broglie wavelength of the particle is much larger 
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FIG. 7. (Color online) Zoom-in of the phase diagram (see 
Fig. [2]) for blue (circles) and red (squares) speckles at the 
bottom of the band. The shape is qualitatively similar to 
the one obtained for spatially correlated speckles in Ref. 0, 
pointing out the crucial role played by the on-site potential 
distribution P{V). 

than the correlation length of the potential, and it is 
reasonable to expect speckles to essentially behave like 
(5-correlated potentials. 

In Fig. [7] we reproduce the portion of the phase dia¬ 
gram near E = —6, calculated in Sec.HIlfor blue and red 
speckle potentials. Several features of previous numeri¬ 
cal calculations [3] for speckle potentials with isotropic 
correlations are indeed recovered: 

• At very low Vq -C E^, for both blue and red- 
detuned speckles, the mobility edge lies below the 
average potential. 

• For larger Vq and blue-detuned speckle, the mobil¬ 
ity edge goes to lower and lower energy. 

• In contrast, for a red-detuned speckle, the mobility 


edge has a turning point and becomes larger than 
the average potential energy. 

However, we stress that a strict mapping of the An¬ 
derson model on the behavior of a massive particle is not 
possible for realistic speckle potentials. Indeed, as dis¬ 
cussed in Ref. 0, the correlation function of the speckle 
potential has a long-range tail (decreasing not faster than 
1/r^ at large distance). As a result, the integral of the 
correlation function diverges, implying that the “white 
noise” limit, which could be associated to a purely un¬ 
correlated potential, does not exist. Thus, even if the 
qualitative behaviors are similar, quantitative differences 
are expected. 

VI. SUMMARY AND CONCLUSION 

In conclusion, we have carried out a thorough analy¬ 
sis of the 3D Anderson model for uncorrelated random 
potentials obeying the Rayleigh probability distribution 
Eq. dl}. We have shown that the asymmetry in the mo¬ 
bility edge of cold atoms exposed to blue and red speck¬ 
les, as recently found numerically in spatially-correlated 
isotropic potentials [3], is directly related to the asym¬ 
metry of the Rayleigh distribution under the change 
V —>■ —V. Finally, our work points out the crucial impor¬ 
tance of a correct evaluation of the single-particle Green’s 
function in speckle potentials, as a necessary condition to 
apply the self-consistent theory of localization. We hope 
that our results will stimulate further work to generalize 
the coherent potential approximation to spatially corre¬ 
lated speckle potentials. 
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